< 

(N 



7~okyo A^etropolitan University 

Nuclear Theory Group, Department of Physics 



TMU-NT940803 

^ ■ |iep-ph/94083gg 



August 1994 



Heavy Quarkonium in 
the dual Ginzburg-Landau theory of QCD 



> 
oo 
cn 

^ ■ Katsuhiko SuzukiQ and Hiroshi TokiQ 



Department of Physics, Tokyo Metropolitan University 
Hachiohji, Tokyo 192-03, Japan 



ix. 

> 
X 

' Abstract 

We study the masses and the leptonic decay widths of heavy quarkonium in the dual 
Ginzburg-Landau (DGL) theory of QCD, in which the abelian monopole condensation 
plays an essential role for color confinement. The effect of color screening due to the light 
quark pair creation is introduced by the infrared momentum cutoff in the gluon propaga- 
tor. We find that the color screening effect is important to reproduce the properties of the 
heavy quarkonium, in particular the leptonic decay width. We also discuss the underlying 
systematics of the heavy quark spectroscopy in the DGL theory. 
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1 Introduction 

It is our strong desire to understand the mechanism of color confinement and to 
construct a low energy effective theory of QCD. Without the knowledge of confinement, we 
can not describe hadrons and nuclei in terms of the fundamental quark degrees of freedom 
as well as can not predict the quark-gluon phase transition under extreme conditions. The 
abelian projection of QCD proposed by 't Hooft is one of the most promising method to 
construct an effective theory of QCD with quark confinement . In this approach, QCD 
reduces to an abelian [f/(l)]^ = U{1)1 x f/(l)g theory with the gauge field singularity, 
which is identified as the magnetic monopole. Recently, this idea is formulated by the 
Kanazawa group as a dual Ginzburg-Landau (DGL) theory of QCD[|2|. By using the DGL 
theory, the quark- ant iquark potential is derived within the static source approximation]^, 
^, which may be realized in the heavy quarkonium. The resulting potential shows the 
linear quark confinement when the QCD monopole condensation takes place. In addition, 
lattice QCD simulations in the maximal abelian gauge|Q indicate the condensation of the 
QCD monopole in the confinement phase^, and the abelian dominance for physical 
quantities related with the confinement such as a string tension]^. These results support 
the 't Hooft's conjecture and the validity of the DGL theory of QCD. 

Beyond the quenched (static quark) approximation, roles of dynamical light quarks are 
important. For the description of the light hadron structure, in which chiral symmetry 
and its spontaneous breaking are essential, we must solve the Schwinger-Dyson (SD) 
equation and the Bethe-Salpeter equation with the [f/(l)]^ field equations together. Even 
if for the heavy quark system, the light quark degrees of freedom are crucial, since they 
produce color screening of the confinement force. Due to the spontaneous pair creation of 
light quarks in the fiux tube, the confinement force between a heavy quark and antiquark 
pair is screened. As a result, heavy-light quark mesons are created, and the heavy quark 
potential should be fiattened at the large distance longer than a typical hadronic size. In 
fact, such a color screening is observed in the lattice QCD calculations with the dynamical 
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fermion, where the confinement potential shows a saturation for r > 1 fni[p. Recently, the 
color screening effect is incorporated within the DGL theory by introducing an infrared 
momentum cutoff in the non-perturbative gluon propagator 0. Here, we study how the 
color screening affects the physics of the heavy quarkonium using the DGL theory of 
QCD. 

The spectroscopy of heavy quarkonia, i. e. cc and bb, is extensively studied within the 
quark potential models0. Such a phenomenological potential assumes the short range 
coulomb and the long range linear confinement force including the the color spin-spin, 
spin-orbit, and tensor interactions, and gives successful descriptions of heavy quarkonia. 
However, those studies reveal the following difficulties; 

1) The leptonic decay widths of cc S-wave excited states, '?/'(3685), ?/'(4040), ■ • •, are 
systematically overestimated in any potential model[P, |lOl, though those models reproduce 
the width of the ground state J/ip as well as the bottomnium leptonic decay width. (See 
Table 2 in ref. 0.) 

2) Conventionally, ^/'(4160) is assigned to the ^Di state of the charm quark- ant iquark 
bound state 0, and therefore the e"'"e~ annihilation width of this state is expected to be 
small, in comparison with the width of S'-state. However, the experimental value is quite 
large, and is the same order of magnitude as those of the S'-states. 

In this paper, we demonstrate how these difficulties are overcome by the introduction 
of the color screening in the DGL theory. For simplicity, we ignore the spin-dependent 
interaction, and consider the spin averaged mass spectrum. In order to include the spin 
and the angular momentum dependence, we must introduce the additional spin-dependent 



forces within the DGL theory as discussed in ref. [11 



2 Heavy quark potential with color screening in DGL 

We briefiy show the heavy quark potential in the DGL theory, and the introduction 
of the color screening effect to the DGL potential. In the framework of the DGL theory. 
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the infrared effective lagrangian of QCD is given by0, 



Cdgl = ^P[tl^{d, + ieA^-H)~m]^-^[n-{dAA)]^-^[n-{dAB)f 
^ [n-{dA A)]\n ■* {d A B)]x + 7^[n ■ {d A B)]\n ■* {d A A)]x 



+ T.Wd,-ge^-BM'-\i\Xa\'~v'f] , (1) 



a=l 



where ip is the quark field with the current mass m, A^ = (Aj^, A^) the abehan gluon 
field with the color operator H = (A'^/2, A^/2), and -B^ = {B^,B^) the corresponding 
dual field. The monopole field Xa has a magnetic charge gSa with the Dirac condition 
eg = 477, and is arbitrary constant vector. The QCD monopole self-interaction is 
phenomenologically introduced in the last term of (|I]), and parameters, A and v, are 
determined later. The DGL lagrangian (|l|) has [f/(l)^]^ x [[/(l)™']^ symmetry, but and 
B^ are not independent. If the magnetic monopole condensation takes place in the QCD 
vacuum, i.e. the magnetic [f/(l)™]^ symmetry is spontaneously broken, the monopole and 
the dual gluon fields acquire their masses, m-^ = 2\/Xv and = y/Sgv. As a result, the 
color electric charge is confined by the dual Meissner effect. Integrating (|ID over 5^, one 
gets a diagonal (abelian) gluon propagator |^, 



D.Ak) = -^{^.. + (ae - 1)^} + (;^^.^(^) ' (2) 

where 

Here, is the dual gluon mass due to the abelian monopole condensation. In the static 
(heavy) quark-antiquark system, we assume the vortex solution of the type-II supercon- 
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ductor. Thus, (0) leads to the potential as0, 



^i'') = ^Yukawa('^) + ^Linear (0 (3) 

q2 Q-niBr 

^Yukawa(^) = "4^^— 

^Linear(^) = H ^ ) ^, (5) 

where is the sum of the diagonal color charge, = Q3 + Ql = If the magnetic 

f/(l)^ X f/(l)™ symmetry were not broken, i.e. dual gluon mass = 0, the potential 
would be a pure coulomb type. The empirical value of the string tension ~lGeV/fm is 
reproduced with ~ O.SGeV and ~ IGeV. 

Up to now, we work with the quenched approximation, and do not deal with an effect 
of the color screening. In the presence of the dynamical light quarks, the string of the 
hadron flux becomes unstable against the spontaneous qq pair creation. This means that 
the flux tube is divided into two pieces when the distance between the valence quark and 
the antiquark becomes larger than a certain critical value Rsc ~ Ifm- Indeed, the results 
of the lattice QCD simulation with the dynamical light quarks show that the qq potential 
becomes flattened for r > lfmp|. 

In the DGL theory, we introduce the color screening effect by the following simple 
replacement in the gluon propagator, as done in ref. 0; 



Here, the parameter e plays an infrared momentum cutoff in the gluon propagator, and 
removes the double pole structure j^^^^, which plays an essential role for the long range 
quark correlation. In principle, this infrared cutoff e can be obtained by calculating the 
quark polarization diagram after getting the non-perturbative quark propagator, though 
it is hard to handle it. In ref. [^, the energy of the created qq pair was estimated 
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< 2Eq > ~ 850 MeV by using the Schwinger mechanism. Since the energy of created pair 
is supphed by cutting off the hadronic string, the critical distance Rsc niay be represented 
as k Rsc ~ < ^Eg >, where k is the string tension ~ IGeV/fm. Then, we may take the 
critical distance of order of Rsc ~ 1 /e ~ Ifm. In the present study, we treat the infrared 
cut e as a free parameter of the model; e ~ lOOMeV. Using (j^), we get the modified heavy 
quark potential with the color screening effect. 



\SC 
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The resulting quark potential for e=0, 100, and 200 MeV is shown in Fig.l. The short 
range Yukawa part is not modified by the color screening. Note that, in the long distance 
limit; r ^ 1/e, the quark potential tends to be constant, and shows a saturation behavior 
which agrees with the lattice QCD calculation. 



^ 0" T inear ' H 2 2 = constant 

^i^mear g^^ ^ ^2 _ ^2 



As shown above, the screening of the confinement force is naturally incorporated in the 
DGL theory of QCD. We remark here that the role of the infrared cutoff e is crucial to solve 
the SD equation for the light quarks, since it removes the infrared double pole structure 
in the gluon propagator]^. We next consider the heavy quarkonium spectroscopy in the 
DGL theory, and study the effect of the color screening on the heavy quarkonium. 

Our present study shares some similarities with the recent work of Yu-bing et a/. ||12||, 
where the color screening effect for the charmonium is studied in terms of the error 
function type confinement force. However, our treatment is fully based on the DGL 
theory, which has a clear link with QCD by virtue of the abelian projection procedure. 
The infrared cutoff parameter e is not merely a phenomenological one, but calculable in 
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the DGL theory. Moreover, we can also study the hght hadron physics in terms of the 
gluon propagator (|^), which enables us to construct a consistent description of light and 
heavy quark systems. In fact, we will show later that the infrared cut e, which is fixed 
by the masses and the leptonic decay widths of the heavy quarkonium, also reproduces 
correct values of light hadron properties such as the pion decay constant and the quark 
condensate < qq >. 

3 Mass spectrum and leptonic decay width 

We calculate the mass spectra of the heavy quarkonium in terms of the quark potential 
obtained in the previous section. In order to get the energy spectrum of bound states, 
we use the diagonalization procedure by using the harmonic oscillator wave function. We 
neglect the spin dependence of the spectrum for simplicity, and discuss the spin averaged 
spectrum, since the energy shift due to the spin dependent interaction is of the order of 10 
MeV due to the large quark mass. First, we show in Fig. 2 the color screening dependence 
of the charmonium spectrum. Here, we use the parameters itlb = 500 MeV, = 1458 
MeV, e = 5.45, the charm quark mass 1289 MeV, and the infrared cutoff e = 0, 100, 
and 200 MeV. It is clear that the IS* ground state (corresponds to J/ip and rjc) is almost 
independent of the choice of e, whereas the spectrum of S-wave excited states and P-, 
D-states receive substantial contributions. Since the change of the confinement potential 
due to the color screening is appreciable for r > Ifm, the IS" state is not affected due 
to its small radius. The 2S, IP, and ID become loosely bound by the reduction of the 
confinement force, and then the energies of these states decrease. To see this behavior 
clearly, we plot in Fig.3 the wave functions of IS", 25", IP, and ID states with e = 0, 100, 
and 200 MeV. Apparently, the IS* wave function is unchanged as expected. Other states 
show considerable e-dependence, in particular, the value at the origin |\I'(0)| decreases as e 
increases. The similar behavior is also found in Fig.4, where the e-dependence of the root 
mean square radii and |^"^|°j|2 are shown. The root mean square radii are normalized 
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to their e = values. The radii of those states at e = MeV are < >\^g = 0.46 fm, 
< >Jp = 0.72 fm, < >2^g = 0.90 fm, and < >J{^ = 0.93 fm. The color screening 
effect makes the sizes of bound states large, and decreases the values at the origin of the 
S'-wave excited states. This behavior is of great significance for the discussion of the e^e~ 
annihilation widths. 

In order to study the color screening effect in detail, we compare two cases; e = 
and e = 100 MeV. For the model parameters, we fix the masses of the dual gluon 
and the magnetic monopole, the gauge coupling e, and the charm quark mass to re- 
producing the masses of IS", IP, 25*, and ID of the charmonium; mB=500(516)MeV, 
m^=1524(1704)MeV, e=4.8(5.2), and m^= 1278(1253)MeV for e=0(100)MeV. In both 
cases, the mass spectrum is well reproduced, as tabulated in Table 1. Note that '?/'(4160) 
is assigned to 4S'-state in the e = 100 MeV case, which is different form the conventional 
2D assignment. The masses of S'-wave excited states decrease due to the screening of the 
linear confinement potential so that above assignment becomes possible. 

The e^e~ annihilation decay width is given by the following formula including the 



QCD correction |T3 



r„.(r^.-.-) = l^|*„.(o)P(i-15^) . (9) 



To avoid the ambiguity of the choice of the QCD correction, we concentrate on the ratio 
of nS to IS' decay width. This ratio corresponds to 



From the comparison with experiment, most of the theoretical calculations overestimates 
the decay widths for 3S, 4S states Our model calculations with e = MeV gives also 
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larger values than the experimental data as found in the third column of Table 2. However, 
in the case e = 100 MeV, the decay width ratio is fairly reproduced as shown in the fourth 
column, since the color screening effect reduces the values of |^E'(0)| as already depicted in 
Fig. 3. Particularly, the results of two cases for ip{4:160) show a striking difference. In the 
case of e = 100 MeV, the calculated width is in a good agreement with data due to the 
45'-state assignment of ^/'(4160). On the other hand, the leptonic decay width vanishes 
in the e = MeV case, since the value of the 2D wave function at the origin is zero. Of 
course, inclusion of the tensor interaction in the potential causes the S-D mixing of the 
wave function, and hence the width of the e = MeV case becomes non-zero. However, 
such a S-D mixing effect gives a still much smaller width than experimental data (about 



10 % at most) within the realistic tensor force [PD 



It is quite interesting to note that the use of e ~100 MeV also leads good results 
for light quark properties such as the pion decay constant /^r ~ 80 MeV and the quark 
condensate < qq >b.gi^ (— 200MeV)^, which are obtained by solving the Schwinger- 
Dyson equation|P, 15|. Both light and heavy quark properties are well reproduced in the 
DGL theory with the color screening. This fact indicates a close connection between the 
confinement mechanism and the chiral dynamics of the light quark physics. 

We finally address a systematics of the heavy quark spectroscopy using several type 
of the potential. In Fig.5, we show the charmonium mass spectra in terms of the linear, 
the linear plus coulomb, and the DGL potential with the experimental data. The linear 
potential provides already quite a good description of the experimental spectrum. The 
coulomb attraction is, however, needed to bring down especially the IS" state relative to IP 
state. The higher 5* states are overestimated with the linear plus coulomb potential. We 
note that tp{4160) was identified to the 2D state due to this agreement with calculation |p. 
If we adopt now the DGL potential with the color screening, the higher S states could be 
brought down largely. This fact makes another assignment of '?/'(4160) to the AS state ||16||. 
If this is the case, we can solve the puzzle of the large leptonic decay width of iplAlQO) as 
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already done in Table 2. 

Here, we comment on the properties of bb in the DGL theory. We get reasonable 
agreements with experimental data using the bottom quark mass = 4734 MeV for 
e—0 MeV. The calculated root mean square radii are < >]^g — 0.20 fm and < >\p 
— 0.41 fm for e = MeV, which are much smaller than those of the charmonium. Since 
the DGL theory is designed as the infrared cfTcctivc theory of QCD, it is not adequate to 
discuss such a short range structure < :>^/'^r^ 0.2 fm ~ lGeV~^. Glearly, the effect of 
the color screening is not so important in the bottomnium case due to those small radii. 
Therefore, we do not discuss the bottomnium spectroscopy here in terms of the DGL 
theory. In order to investigate the bb, we must consider the short range gluon exchange, 
i.e. coulomb type interaction, in which the role of non-diagonal gluon is crucial. We have 
to develop a theory, which connects smoothly the long distance physics discussed here 
and the short range perturbative physics. Such a study is beyond the scope of the present 
work. 



4 Summary and discussions 

We have studied the heavy quarkonium properties in the dual Ginzburg-Landau theory 
of QGD, which generates the color confinement due to the abelian monopole condensation. 
The color screening of the confinement potential due to the light quark pair creation 
has been incorporated by introducing an infrared momentum cutoff in the abelian gluon 
propagator. Such an infrared momentum cut controls the quark long range correlation. 

We have found that the DGL theory gives a good description of heavy quarkonia with 
e ~100 MeV. The screening of the confinement force reduces the leptonic decay width of 
the nS{n > 2)-states. This potential brings down also the higher S states, and make it 
possible to identify ■?/'(4160) as the 45* state instead of 2D state. This identification solve 
the problem of the large e+e" decay width of ■^(4160). 

Further studies on the DD decay of the charmonium and BB decay of the bottomnium 
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are of special interest within the DGL theory. For those processes, the quark pair creation 
model|TU, |T7| can reproduce the experimental data. In such a study, however, a probability 
of the light quark pair creation from the vacuum is treated as a free parameter. In the DGL 
theory, the rate of the light quark pair creation can be estimated by using the Schwinger 
mechanism!]^. Realistic calculations of the pair creation process are under investigation. 

It is very interesting and important to note that the choice of e ~100 MeV also 
reproduces the pion decay constant and the quark condensate in the light quark physicsQ. 
Commonly, the low energy hadron physics is believed to be dominated by the chiral 
dynamics, whereas the heavy quark system is governed by the linear confinement potential 
with the perturbative coulomb part, in which the chiral symmetry is less important. 
Hence, we have no clear connection between the confinement and the physics of the chiral 
symmetry. On the other hand, the lattice QCD at finite temperature tells us that the 
deconfinement transition and the chiral transition take place at the same temperature, and 
suggests that the quark confinement and the chiral dynamics should be strongly correlated. 
In the DGL theory of QCD, the quarks are confined by the dual Meissner effect due to the 
abelian monopole condensation. Here, we have shown that the heavy quark spectroscopy 
as well as the light hadron properties are reproduced by the same fundamental parameter 
of the DGL, i.e. the monopole mass, the dual gluon mass, the gauge coupling, and 
the infrared momentum cut e. Such agreements support the DGL theory as the effective 
theory of QCD, which provides an unified treatment of the confinement and the dynamical 
chiral symmetry breaking. These findings motivate us to do further investigations of 
various low energy phenomena in terms of the DGL theory of QCD. 



10 



References 

[1] G. 't Hooft, Nucl. Phys. B190 (1981) 455 

[2] T. Suzuki, Prog. Theor. Phys. 80 (1988) 929; 80 (1989) 752 
S. Maedan and T. Suzuki, Prog. Theor. Phys. 81 (1989) 229 

[3] H. Suganuma, S. Sasaki and H. Toki, RIKEN AF-NP-164 (1994), |hep-ph/9312350| 
H. Toki, H. Suganuma and S. Sasaki, in Proceedings of International Symposium on 
Spin-Isospin Responses and Weak Processes in Hadrons and Nuclei, eds. T. Suzuki 
et at, to appear in Nucl. Phys. A (1994) 

[4] For a review, T. Suzuki, Nucl. Phys. B (Proc. Suppl.) 30 (1993) 176, and references 
therein 

[5] A.S. Kronfeld et al, Phys. Lett. B198 (1987) 516 

A.S. Kronfeld, G. Schierholz and U.-J. Weise, Nucl. Phys. B293 (1987) 261 

[6] S. Hioki et al, Phys. Lett. B285 (1992) 343 

[7] S. Hioki et al, Phys. Lett. B272 (1991) 326 

[8] E. Laermann et al, Phys. Lett. B173 (1986) 437 
R. Gupta et al, Phys. Rev. D44 (1991) 3272 
W. Sakuler et al, Phys. Lett. B276 (1992) 155 

[9] E. Eichten et al, Phys. Rev. D21 (1980) 203 

S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189 

For a review, Quarkonia, ed. W. Buchmiiller, North-Holland 

For a comparison of models, see D. Besson and T. Skwarnicki, preprint CLNS 93/1239 
(1994), and references therein 



11 



[10] K. Heikkila, N.A. Tornvist, S. Ono, Phys. Rev. D29 (1984) 110 

[11] M. Baker, J.S. Ball and F. Zachariasen, Phys. Rev. 44 (1991) 3949,i&idD47 (1993) 
3021; Phys. Lett. B283 (1992) 360 

For a review, M. Baker, J.S. Ball and F. Zachariasen, preprint CALT 68-1887 (1994) 

[12] Dong Yu-bing et al, Phys. Rev. D49 (1994) 1642 

[13] R. Van Royen and V. Weisskopf, Nuovo Cimento A50 (1967) 610 
R. Barbieri et al, Phys. Lett. B57 (1975) 455 

[14] Particle Data Group, K. Hikasa et al, Phys. Rev. D45 (1992) SI 

[15] S. Sasaki, H. Suganuma and H. Toki, TMU preprint (1994) 

[16] C. Quigg and J. Rosner, Phys. Rep. 56 (1979) 167 

[17] A. Lc Yaouanc et al, Hadron Transitions in the Quark model; Gordon and Breach 
(1988), and references therein 



12 



Table Captions 

Table 1 : The charmonium mass spectrum in the DGL theory with color screening. 



The experimental data|Tj| with the corresponding states are shown in the first column. 
Theoretical results with e = MeV and 100 MeV are depicted in the second and the third 
columns, respectively. Note that we do not consider the spin-dependent interaction, and 
calculated values are the spin averaged ones. 



Table 2 : The ratio of nS to 15* leptonic decay width of the charmonium in the DGL 
theory. The experimental data[|l^ are shown in the second column. The results with e = 
MeV and 100 MeV are depicted in the third and the fourth columns, respectively. 
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Table 1 



State 


e = 


e = 100 




1,5(3074) 


15(3066) 


Xcj{^ 3500) 


1P(3492) 


1P(3507) 


^(3685) 


25(3664) 


25(3644) 


^(3770) 


1^(3781) 


1D(3775) 


^(4040) 


35(4108) 


35(3999) 


^(4160) 


2D(4190) 


45(4256) 


V'(4415) 


45(4492) 


55(4469) 



Table 2 



State 


Exp 


e = 


e = lOOMeV 


rv)(3685)/rj/^ 


0.40 ± 0.04 


0.48 


0.38 


r^(404o)/rj/^ 


0.14 ±0.03 


0.34 


0.22 


rv)(4i6o)/rj/v, 


0.14 ±0.04 


0.0" 


0.18'' 


rv,(44i5)/rj/^ 


0.09 ±0.02 


0.27 


0.08 



a) In this case ^/'(4160) is assigned to the 2D state due to the mass spectrum. 

b) In this case ■0(4160) is assigned to the 45 state due to the mass spectrum. 
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Figure Captions 

Fig.l : The heavy quark potential in the dual Ginzburg-Landau theory with the color 
screening. The solid curve represents the case of e = MeV (no screening). The dash- 
dotted and the dashed curves denote the e = 100 and 200 MeV cases, respectively. 

Fig. 2 : The infrared cutoff e dependence of the charmonium spectrum. 

Fig.3 : The wave functions of the cc bound states for 15", 25", IP, and ID states calculated 
in the DGL theory with the color screening. The results with e = 0, 100, and 200 MeV 
are shown by solid, dash-dotted, and dashed curves, respectively. 

Fig.4 : The infrared cutoff e dependence of the root mean square radii and the values 
of the wave functions at the origin. The root mean square radii normalized to the e—0 
MeV value, are shown by the dashed curves for 15", IP, ID, and 25" states. |^"^'|Qjj2 is 
also shown by the sohd curves for the 25" and 35" states. 

Fig.5 : The charmonium mass spectra in various potential models; from left to right the 

linear potential, the linear plus coulomb potential, the DGL theory with color screening, 
and experimental spectrum. In the experimental spectrum, the state at 4160MeV is either 
assigned to the 2D state (a) and 45" state (b). 
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Fig.l 
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